Density Functional Study of Excess Fe in Fei+^^Te: Magnetism and Doping 
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The electronic and magnetic properties ol the excess Fe in iron telluride Fe(i_|_j.)Te are investigated 
by density functional calculations. We find that the excess Fe occurs with valence near Fe"*", and 
therefore provides electron doping with approximately one carrier per Fe, and furthermore that the 
excess Fe is strongly magnetic. Thus it will provide local moments that interact with the plane 
Fe magnetism, and these are expected to persist in phases where the magnetism of the planes 
is destroyed for example by pressure or doping. These results are discussed in the context of 
superconductivity. 

PACS numbers: 74.25. Jb,74.25.Kc,74.70.Dd 



Recently, iron chalcogenides a-FeSe and a-FeTe, an- 
other new family of Fe-based superconductors have been 
reported.— i^i^i^i^i^ii The superconducting transition tem- 
perature Tc has increased from initial 8 KA to 14S. (15.3^) 
K with appropriate Te substitution, and 27 K at high 
pressures (1.48 GPa).^ While the presently known max- 
imum critical temperatures are lower than in the Fc-As 
families (^iSiiiSiiiii^ this binary system has drawn consid- 
erable attention due to the apparent simplicity of the 
structure, the fact that it is As free, and the fact that 
large crystals of Fei+j;(Se,Te) can be grown enabling de- 
tailed characterization by neutron and other measure- 
ments. These compounds the a-PbO structure, which 
consists of a c-axis stack of FeTe sheets, with each sheet 
consisting a square planar layer of Fe, tetrahedrally co- 
ordinated by Te, similar to the FeAs sheets of LaFeAsO 
or LiFeAs. In fact from a structural point of view these 
compounds are very similar to LiFeAs, with As replaced 
by a chalcogen and the Li replaced by a site with a low 
partial filling of excess Fe. According to literature, these 
compounds always form with excess Fej -'^^'^'^i^^'^^i^^'^^ 

Electronic structure calculations for the stoichiomet- 
ric iron chalcogenides, FeATji^ show electronic structures 
and Fermi surface topologies very similar to those of the 
other Fe-based superconductors. ^^'^-^i^^'^'^i^^'^^i^^ There is 
a general proximity to magnetism, especially in FeTe, as 
well as a substantially nested Fermi surface, which fa- 
vors a spin density wave (SDW) instability at the 2D 
(tt,tt) point. While the mechanism for superconductiv- 
ity in the Fe-based superconductors is yet to be estab- 
lished, there is a strong association between the occur- 
rence of the SDW and superconductivity in the phase dia- 
grams, with superconductivity generally occurring when 
the SDW is destroyed either by doping or by pressure. 
The SDW is observed in most of the undoped Fe-As su- 
perconducting materials, and is accompanied by a lattice 
distortioni22i22ii^i^i^ For the chalcogenides, a structural 
distortion with decreasing temperature was detected in 
FeSe(i_2;|^ and superconductivity was found to be close 
to magnetic instability in Fe(Se(i_2;)Tej;)o.82 (the for- 
mula does imply chalcogen vacancies but reflects excess 
Fe).® Furthermore, Fei+a;Te is reported as magnetic, with 
properties depending on stoichiometry in several older 




FIG. 1: (Color online) Structure (2x2 supercell of the tetrag- 
onal a- FeTe with one excess Fe) used to simulate Fei.i25Te 
(composition FegTeg). The iron in Fe-Te layers are denoted 
Fei and excess iron as Fe2. 



papers. Bao et alM^ based on neutron results suggested 
a more complex incommensurate antiferromagnetic order 
for the Fe(Se(i_a;)Tea;) system than in the Fe-As based 
SDW phases. On the theoretical side, magnetism driven 
by Se vacanciesr^ non-coUinear— and bi-coUinear— an- 
tiferromagnetic state have been suggested. 

Here we report supercell calculations investigating the 
role of the excess Fe focusing on Fei+ccTe. We find that 
as might be expected, excess Fe donates charge to the 
FeTe layers, acting as an electron dopant. Interestingly, 
it occurs with a valence near Fe"*" with each Fe donating 
one carrier. Furthermore, there is a very strong tendency 
towards moment formation on the excess Fe. These mo- 
ments will then interact with the magnetism of the FeTe 
layers, perhaps complicating the magnetic order. They 
would also be expected to persist into the regime where 
FeTe magnetism is suppressed by doping or pressure, per- 
haps extending the range of magnetic order in the phase 
diagram, and providing pair breaking in the supercon- 
ducting state. 
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FIG. 2; (Color online) Calculated electronic total and partial 
DOS for non-polarized Fei.i25Te. 



The electronic structure and magnetism calculations 
were performed with the projector augmented wave 
method^ as implemented in VASP code.— The gen- 
eralized gradient approximation^ was employed for the 
exchange-correlation functional. A kinetic energy cut- 
off of 268 eV and augmentation charge cutoff of 511 eV 
were used to obtain converged energy (within 1 meV). 
To simulate the partially occupied excess Fe, we used a 
2x2 supercell of a-FeTe (two formula per cell) with one 
Fe atom (labeled as Fe2) placed at the 2c (G.5,G,z) site, 
as shown in Fig. [TJ This corresponds to a stoichiometry 
of Fei.i25Te. The experimental lattice parameters a — 
3.8245 and c = 6.2818 for Fei.i25Te Ail were used in our 
calculations. An 8x8x10 grid was used for the fc-point 
sampling of the Brillouin zone , and a denser 16x16x20 
fc-mesh was used for density of state (DOS) calculations. 
The internal coordinates were relaxed to minimize the 
forces to below 0.01 eV/A. The calculated coordinate of 
Fe2 is zpe — 0.703. This is in reasonable agreement with 
the experimental results in Ref. [ij {zpe = 0.692) and 
Ref. i lSl {zF e = 0.721), but significantly higher than that 
in RefTlB [zpe = 0.561). 

We begin by showing that the excess Fe atom is 
strongly magnetic. Fig. [2] shows the DOS for Fei.i25Te 
obtained in a nonmagnetic calculation. The electronic 
states near the Fermi level {Ep) are mostly of char- 
acter of the Fei layers with small contribution from the 
excess Fe2 atom. The result that the Fermi level lies ex- 
actly at a sharp peak of the Fe2 3c? DOS indicates the 
magnetic instability. The calculated Fe2 partial DOS at 



TABLE I: Calculated energy difference (in meV/Fei.i25Te, 
relative to the non-magnetic state for all of Fe) between dif- 
ferent types of magnetic arrangements for Fei in layers and 
excess Fe2. 

Fe2 (non-magnetic) Fc2 (magnetic) 
Fei (non-magnetic) -48.9 

Fei (ferromagnetic) -104.9° -72.0 

Fci (checkerboard antiferromagnetic) -93.1 -151.0 

Fei(SDW antiferromagnetic) -131.0 (-199.6'') -186.7 (-256.6'') 

"Actually, this type of arrangement always converges to an "anti- 
ferromagnetic" order with Fei and Fe2 having opposite spin direc- 
tions. 

''The calculated value with relaxed structure when fully consider- 
ing magnetic orders. 



Ep is 6.2 states/eV/Fe (both spins). Within the Stoner 
theory the magnetism occurs when N{Ef)I > 1, where 
N{Ep) is the DOS at the Fermi level per atom per spin 
and / is Stoner parameter, typically in the range of 0.7 ~ 
0.9 eV for Fe. The large Fe2 3d DOS at £'f easily satisfies 
the Stoner criterion for the magnetic instability. Indeed, 
when considering the spin polarization for the Fe2 atom 
(not Fei layers), the total energy is reduced by 48.9 meV 
as shown in Table |T1 A pseudogap is opened with the 
Fermi level falling into its bottom, as shown in Fig. [31 
The calculated Fe2 DOS at the Fermi level is reduced to 
0.8 states/eV/Fe. The calculated magnetic moment for 
Fe2 is 2.5 fiB- 

Based on integration of the partial Fe2 DOS up to the 
Fermi level and normalization with the total Fe2 DOS 
we find 4.7 electrons in the majority spin states and 2.2 
electrons in the minority spin states. Thus, the excess Fe 
occurs as Fe+ and each excess Fe atom has donates ap- 
proximately one electron to the Fei layer. It may noted 
that Fe+ is a somewhat unusual valence state for stable 
Fe compounds. Here this state is stabilized because of 
a balance between the two Fe sites. Specifically, in stoi- 
chiometric FeTe, Fe is already di-valent, and so the more 
rapid electron doping that would result if the excess Fe 
were di-valent would lead to a more rapid conversion of 
the plane Fe towards Fe+ . This balance between low va- 
lence states for Fe in the plane and excess positions may 
be responsible for the fact that the structure does not 
form at higher excess Fe concentrations. 

For undoped FeTe, the Fermi level is located somewhat 
below the bottom of the pseudogap as shown in Ref. [H. 
The presence of the excess Fe atoms in Fei.i25Te moves 
the Ep up, reducing the total DOS at Ep. However, de- 
spite the electron doping, the total DOS at Ep remains 
relatively high (1.8 states/eV/Fe), which would still put 
the Fei.i25Te close to magnetic instabilities. Our calcu- 
lations show that stripe antiferromagnetic ordering (the 
SDW type) is most stable compared to the nonmagnetic, 
ferromagnetic, and the checkerboard antiferromagnetic 
phases (see Table H]) assuming the fixed structure for the 
nonmagnetic phase in all these calculations. Relaxing 
the structure for the SDW antiferromagnetic state fur- 
ther lowers the total energy by nearly 70 meV. 
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FIG. 3: (Color online) Calculated electronic DOS for 
Fei.i25Te with the moment formation on the excess Fe2 (non- 
magnetic order for Fei layers). 



The magnetic moment of the Fe2 is calculated to be 
2.4 /is, much higher than that for the Fei layers (1.6- 
1.8 uB). The excess Fe's strong magnetism is supported 
by the recent neutron scattering experiment.^* The large 
local magnetic moment of the excess Fe is expected to 
persist even if the SDW antiferromagnetic ordering of the 
Fe layers is suppressed by the doping or pressure, thus 
causing pair breaking in the superconducting phase. 

As we have found in our early work, FeTe has the 
strongest SDW state in the iron chalcogenide familyJ^ 
This is consistent with experimental observations that 
the SDW state is maintained up to high excess Fe con- 
tents, which as discussed correspond to high doping levels 
(with X up to 0.125 in the Fei+2;Te systems). "^'^^ It should 
also be noted that the heavy doping corresponds to a 



large size mismatch between the approximately cylindri- 
cal electron and hole Fermi surfaces. The driving force 
for an itinerant spin density wave is Fermi surface nest- 
ing. As notedj^^ for a large size mismatch the structure in 
the susceptibility around (1/2,1/2) will develop a dip at 
the center, with the maximum therefore moving off cen- 
ter. If the SDW stays stable the ordering vector will then 
become incommensurate. This apparently is the case in 
Fei+a;Te, and may explain the incommensurate SDW ob- 
served in neutron scattering.-'^* We note than an itinerant 
SDW can arise simply within this itinerant framework, 
but would require a complex frustration within a local 
moment picture. 

Turning to the trends, as noted FeTe has a 
stronger tendency towards magnetism than FeSe and 
the arsenides, but also still shows signatures of spin 
fluctuationsfi^ As such, within a scenario where super- 
conductivity arises from pairing due to spin fluctuations 
associated with the Fermi surface nesting, FeTe may have 
particularly high temperature superconductivity if the 
SDW can be suppressed. However, the SDW persists 
up to high doping levelsJ^ A particularly interesting ex- 
periment would be then to destroy the SDW by pres- 
sure and search for superconductivity in the resulting 
paramagnetic phase. Furthermore, the fact that the ex- 
cess Fe in this compound, and presumably the Fei+a:Se 
and Fei+2;(Se,Te) superconductors, has a local moment 
in proximity to the Fe layers offers an interesting oppor- 
tunity for experimental investigation of the interplay be- 
tween superconductivity and presumably pair breaking 
magnetic scattering in the Fe superconductors. 

In summary, we find that the excess Fe in Fci+^^Te 
is strongly magnetic and is also an electron donor, with 
each excess Fe atom donating approximately one electron 
to the Fe layer. 
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